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Numerical simulations by means of Monte Carlo method and finite-size scaling analysis have been 
performed to study the percolation behavior of linear /c-mers (also denoted in the literature as rigid 
rods, needles, sticks) on two-dimensional square lattices L x L with periodic boundary conditions. 
Percolation phenomena are investigated for anisotropic relaxation random sequential adsorption of 
linear /c-mers. Especially, effect of anisotropic placement of the objects on the percolation threshold 
has been investigated. Moreover, the behavior of percolation probability Rl(p) that a lattice of 
size L percolates at concentration p has been studied in details in dependence on k, anisotropy and 
lattice size L. A nonmonotonic size dependence for the percolation threshold has been confirmed 
in isotropic case. We propose a fitting formula for percolation threshold p c = a/k a + frlog 10 k + c, 
where a, 6, c, a are the fitting parameters varying with anisotropy. We predict that for large Ai- 
mers (k 1.2 x 10 4 ) isotropic placed at the lattice, percolation cannot occur even at jamming 
concentration. 

PACS numbers: 64.60.ah, 64.60.De, 68.35.Rh, 61.43.Bn 

Keywords: percolation, disordered systems, Monte Carlo Simulations 



I. INTRODUCTION 

Percolation deals with the properties of disordered me- 
dia. Such media can be composed of the objects placed in 
a space. The objects can connect with each other and 
form clusters. If object concentration is large enough, 
infinitely large cluster occurs. Such a concentration is 
known as a percolation threshold. The properties of me- 
dia are considerably different below and above percola- 
tion threshold. If objects are placed in a space purely ran- 
domly, the percolation is called random or Bernoulli per- 
colation. Moreover, different correlations or constrains 
may be applied to the space distribution of the objects. 
The media composed in such a way may be partially dis- 
ordered and anisotropic. Very often, a discrete space (lat- 
tice) is utilized to simplify consideration. In this case, the 
cluster-forming objects are sites of the lattice. Percola- 
tion of the point objects (singly occupied site) on differ- 
ent lattices in plane and multidimensional space is more 
intensively studied. Percolation of the objects occupying 
several nearest sites is studied much worse. The examples 
of such objects are linear, cyclic and branched /c-mers, 
i.e. k nearest sites. The huge amount of publications is 
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devoted to both theoretical and applied aspects of per- 
colation (see, e.g., [lH3|). During the past few decades, 
percolation of the anisotropic penetrable and impenetra- 
ble objects (rods, sticks, linear fc-mers, ellipsoids etc.) 
has been intensively investigated. In our overview, we 
restrict ourselves to the works devoted to the percolation 
of linear objects on a lattice. 

Mainly, the studies are devoted to the isotropic prob- 
lem on a square lattice when the fc-mers with horizontal 
and vertical orientations are deposited with equal prob- 
ability. A computer-simulation model for linear k-mers 
(k = 1 . . . 20) showed that percolation threshold p c de- 
crease with increasing of the chain- length k as 1/fc 0,5 
The percolation exponents (order parameter, suscepti- 
bility, and correlation length exponents) seem to remain 
unchanged. 

The study of the percolative properties of systems gen- 
erated by a random sequential adsorption (RSA) of k- 
mers (k = 1 . . . 40) have been performed by Leroyer and 
Pommiers [5]. They demonstrated that as the segment 
length grows, the percolation threshold p c decreases, goes 
through a minimum and then increases slowly for large 
k (k > 16). 

Later on, Kondrat and P§kalski [6j extended the stud- 
ies percolation and jamming of the same problem to the 
k-mer length in the interval k = 1 . . . 2000. The authors 
showed that the jamming threshold decreases monotoni- 
cally approaching the asymptotic value of pj = 0.66±0.01 
at large k and percolation threshold p c is a nonmonotonic 
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function of the length /c, showing a minimum for a certain 
length of the k-mers (k = 13). However, these results for 
very large needles cannot be treated as accurate because 
of moderate size of the studied lattices (L < 2500) and 
possibility of large finite-size corrections. 

The details of the monotonic behavior of the perco- 
lation threshold for small k-mer length (k < 15) have 
been widely discussed in the literature |7H10|. Percola- 
tion and jamming phenomena have been investigated for 
k-mer length within the interval k = 1 ... 10 by Van- 
dewalle et al. [7]. The authors conjectured presence of 
some universal connection in the geometry of jamming 
and percolation that resulted in constancy of the ratio of 
percolation and jamming concentration p c /Pj (— 0.62) 
for all sizes of k-mers. The following equation for the 
percolation threshold as a function of k-mer length has 
been proposed 

»- c { 1 -<( L rf)- (1 > 

where C and 7 are the constants. 

Cornette et al. [8] performed the finite-size scaling tests 
and shown that the k-mer problem in all the studied 
cases, belongs to the random percolation universality 
class. They fitted the data for the k-mers (k = 1 . . . 15) 
with following exponential equation: 

p c =p? + nexp(-^, (2) 

where p™ = 0.461 ± 0.001, SI = 0.197 ± 0.02, and k = 
2.775±0.02 are the fitting parameters, p^ is the expected 
value in the limit k — »> 00. 

Recently, these problems have been extended for par- 
tially ordered k-mer (when the particles with horizontal 
and vertical orientations can be deposited with of unequal 
probability) \lM\M- 

The effects dimer alignment on per- 
colation and jamming phenomena on a square lattice has 
been investigated by Cherkasova et al. [llj . The influ- 
ence of dimer alignment on the electrical conductivity 
has been examined, too. The effect of k-mer alignment 
on the jamming threshold has been extensively examined 
for the k in the interval 1 . . . 256 [12]. The percolation 
behavior for the k-mer length in the interval k = 1 ... 15 
has been studied recently by Longone et al. [l3|. Only 
two particular cases have been studied in the work, i.e. 
the isotropic case and the completely ordered case (all 
k-mers are aligned along the given direction). In both 
cases, the percolation threshold is monotonic decreasing 
function of the k-mer length k. 

The numbers of a numerical studies have been recently 
devoted to the anal ysis of equilibrium properties in sys- 
tems of k-mers [14j-ll7|. The equilibrium systems have 
been simulated using the deposition-evaporation dynam- 
ics. The studies showed existence of a orient at ionally 
ordered phase (nematic phase) for long fc-mers. The uni- 
versality class for the percolation and isotropic-nematic 
phase transition have been found to be the same as of 



the random percolation and Ising models. The non- 
monotonic size dependence has been observed for the per- 
colation threshold of unaligned /c-mers, it goes through 
a minimum at k ~ 5, and asymptotically converges to- 
wards a definite value p c ~ 0.54 for large fully aligned 
k-mers [18] . It has been interpreted as a consequence of 
the isotropic-nematic phase transition occurring in the 
system for large values of k. 

Except pure theoretical interest, such considerations 
may have different applications. For instance, the perco- 
lation approach is suitable to describe physical and chem- 
ical properties of monolayers formed during adsorption of 
the polymer chains [19] . Another possible application is 
connected with the nanotechnologies (see, e.g., [20]). Re- 
cently, the current progress on the production of aligned 
single- walled carbon nanotubes (SWCNTs) has been re- 
viewed by Ma et al. [21]. The semiempirical theories 
of composites containing randomly oriented anisotropic 
inclusions (needle, prolate or oblate spheroid, sphere, 
or disk) have been developed and they are useful for 
prediction of effective electrical or thermal conductivi- 
ties of multi- walled carbon nanotube composites [2(1 122I - 
[25(l . The first experiments evidenced the lowering of the 
threshold respective to isotropic systems [26] . The exper- 
iments for random stick patterns obtained by photolitho- 
graphic techniques supported the universality hypothesis 
for 2d systems [27|. The universality concept has been 
also confirmed in experiments with the aluminum film 
containing the insulating ellipsoids with the same direc- 
tion of the major axis |28| . 

This work discusses the percolation behaviour of linear 
k-mers on square lattice with different degree of align- 
ment characterized by order parameter. We try to clear 
the uncertainty in question about the presence or ab- 
sence the nonmonotonic /c-dependence for the percolation 
threshold by studying the systems with k varies from 1 
up to 512. 

In our work, we try to find the answers to the questions 
listed below 

1. Are the Eqs. [1] and [2] valid for the very long linear 
objects or they work only for rather short objects? 

2. How does anisotropic placement of the objects ef- 
fect the percolation threshold? 

The rest of paper is arranged as follow. In Section [III 
we describe our model and the details of simulation. The 
obtained results are discussed in Section HIB We summa- 
rize the results and conclude our paper in Section IIVI 



II. DESCRIPTION OF MODELS AND DETAILS 
OF SIMULATIONS 

The problem of linear k-mers of length 2 n , where n = 
1, 2, . . . , 9, on the square lattices of size L x L has been 
studied. Linear lattice size, L, varies from 100 to 19200 
in different simulations. Periodic boundary conditions in 
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vertical and horizontal directions have been applied, i.e 
percolation on a torus has been considered. 



A. Filling of the lattice by /c-mers 

The relaxation random sequential adsorption (RRSA^ 
model [12] has been used to place the k-mers on a lattice 
In this model, there is an infinitely large reservoir fillec 
with /c-mers oriented with given and fixed anisotropy 
The k-mer is taken from the reservoir and an attemp 
of its deposition is carried out starting from a lattice sit 
selected at random until the object is deposited. In cod 
trast with the conventional random sequential adsorptio: 
(RSA) model, when any unsuccessful attempt is rejectee 
and other object is selected for deposition, RRSA mode 
ensures that anisotropy of the deposit is the same as th 
anisotropy of the objects suspended in the reservoir [12] 

The degree of anisotropy is characterized by the orde 
parameter s defined as 
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where N\ and 7V_ are the numbers of /c-mers oriented ii 
vertical and horizontal directions, respectively. 

For isotropic system, s = 0, the numbers of vertical am 
horizontal /c-mers are the same, and for totally aligne< 
system, s = 1, all k-mers are aligned in vertical directior 
For these two marginal cases, RRSA and RSA models ar 
absolutely identical [l2j |. 

The Mersenne twister random number generator (29) 
with a period of 2 19937 — 1 has been exploited to generate 
positions and orientations of the deposited objects. 



B. Determination of percolation threshold 
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(b) Spirallike wrapping clusters 
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(c) Ringlike wrapping clusters 

Figure 1. (Color online) Percolating clusters of different sorts 
on a plane and on a torus. 



A crossing cluster is determined as a cluster that con 
nects two opposite borders of lattice with open boundary 
conditions. Examples of crossing clusters that percolat< 
along vertical direction or simultaneously along vertica 
and horizontal directions are presented in Figure [TJi. 

A wrapping cluster is determined as a cluster which 
winds (i.e. provides a path of length 2ir) around the 
lattice with the periodic (toroidal) boundary conditions 
along the given direction [30| . The wrapping cluster may 
be either disconnected (spirallike) (Figure [TJ)) or contin- 
uous (ringlike) (Figure Ht), or more complex one. 

From the topological point of view, the spirallike clus- 
ters presented in Figure [TJ) are homotopic to a point, 
i.e, they can be continuously deformed to a point, and 
hence essentially differ from the ringlike clusters shown 
in Figure [T^. From the physical point of view, it is rather 
natural to think that applying periodic boundary condi- 
tions cannot destroy a percolating state existing in plane 
with open boundary conditions. Moreover, it can pro- 
duce a new percolating state because of additional kind 
of symmetry, i.e. translation symmetry. 



In our study, a system is considered as percolating if at 
least one spiral cluster (Figure [Dd) can be found. For def- 
initeness we call it as a problem of physical percolation on 
a torus in contrast with topological percolation when only 
self-connected clusters are treated as wrapping ones |31j . 

The value of threshold concentration may be deter- 
mined by calculation of the probability Rl(p) for a clus- 
ter to cross a square lattice of L x L sites, if the bound- 
ary conditions are open, or to wrap around the peri- 
odic boundary conditions. In the thermodynamical limit 
(L —> oo), this probability is equal to the probability that 
the system percolates (i.e. it tends to the step-function 
and equals below the percolation threshold and 1 above 
it) Hi. 

Since cluster wrapping can be defined in a number of 
different ways (see, e.g., [32]) there are a corresponding 
number of different probabilities Rl' 

1. is the probability of wrapping horizontally 
around the system; 

2. R\ is the probability of wrapping vertically around 
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the system; 

3. is the probability of wrapping around either the 
horizontal or vertical direction, or both; 

4. Rf^ d is the probability of wrapping around both 
directions simultaneously. 

For the square lattices and isotropic problem these 
probabilities satisfy the following relations [32l. l33j: 

R h L = Rl, (4) 

R\ = (Rf + Rf d )/2, (5) 

as well as the inequalities 

Rl nd <R h L < R° L r . (6) 

Relations ([3]), and ([5]) evidence that only two of per- 
colation probabilities are independent. Obviously, for 
an anisotropic system the relation ([!]) cannot hold and, 
hence, there are three independent probabilities. Never- 
theless, for a strong anisotropic systems a spanning or 
wrapping cluster always arises along one direction, say 
vertical, and hence R v = R or , R h = R and . 

The detailed studies have shown d, Q that for the 
specified problem (e.g., for crossing or wrapping clusters) 
and the criterion used, the curves Rl(p) cross each other 
in a unique intersection point i?* located at p = p c in the 
thermodynamical limit (L — >• oo). 

Figure [2] compares R™{p) and Rf^ d (p) dependencies 
for monomer problem (k = 1), i.e. conventional site prob- 
lem and different size of square lattice, L. The results are 
presented for the systems with periodic and open bound- 
ary conditions. 

If the disconnected spiral clusters similar to shown in 
Figure [TJ) are not treated as percolating, the exact ex- 
pressions of R* at percolation threshold, p c , for each of 
the definitions have been deduced [34] from the work 
by Pinson [3lJ. The values of R* presented by New- 
man and Zirf[33,|33 are iT or = 0.690 473 725, iT and = 
0.351642 855. 

The intersection points R* for physical percolation in 
our study are R^ T ~ 0.90 and R* L and ~ 0.98. 

The Rl(p) functions have been estimated by perform- 
ing 1000 independent runs. Percolation concentration 
p c (L) for the lattice of given linear size L filled with k- 
mers at the given concentration p has been determined 
using the fitting function [35] 

R L (p) = (1 + exp (-(p - PcWd))- 1 , (7) 

where a is adjusted constant. 

To extrapolate the estimations of the percolation 
thresholds p c (L) obtained at the lattice of size L to the 
infinite large lattice p c (oo), the usual finite size scaling 
analysis of the percolation behavior has been done. To 
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Figure 2. (Color online) Comparison of R°^ and Rf^ d versus 
p dependencies for monomer problem (k = 1) (physical perco- 
lation) and crossing clusters (with open boundary conditions) 
and different size of square lattice, L. 

perform extrapolation, we used at least three lattices of 
different sizes and scaling relation 

IPc^-PctocOlocL- 1 /", (8) 

where v — 4/3 is the critical exponent of correlation 
length for the 2d random percolation problem [1]. The 
universality of the fc-mers problem has been justified ear- 
lier [8]. In our study, the typical values of lattice size are 
L = 50k, 75/c, lOOfc, 150Jfe, 200/c, 400/c. 

Examples of p c versus L scaling behavior for k = 16, 
s = 0.8 and four criteria (h,v,or,and) are presented in 
Figure 03 

The preliminary studies have shown that in all cases 
the p c (L) scaling is minimal for criterion and and the final 
results on percolation concentration have been obtained 
using the criterium and. To simplify the notation, below 
we omit superscript and where it is possible. 

To avoid very time-consuming computations with the 
lattices of huge size for k = 256 and s = , we used only 
two relatively small lattices L = 50k and L = 75k and 
two different criteria, namely and and or. Intersection 
points (i.e. p c (oo)) extracted from (|5J) for two different 
criteria are almost the same within error bar about 0.001. 

Another special case is L = 512, 8 = 0. Only one 
lattice size L = 37k has been used for rough estimation 
of the percolation threshold. Percolation concentration 
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(a) Isotropic case, k = 2, s = 0.0 (b) Slight anisotropic case, k = 4, s = 0.1 (c) Anisotropic case, k = 16, s = 0.8 

Figure 3. (Color online) Percolation concentration p c versus size of the lattice L for and different criteria for physical percolation. 



has been calculated from the equation R^ d (p c ) = 0.9. 



C. Other details 

Breadth-first search (BFS) algorithm has been applied 
to identify a percolation cluster. BFS seems to be faster 
and more appropriate for the toroidal boundary con- 
ditions than Hoshen-Kopelman (HK76) algorithm (36| . 
The additional tests have shown that results obtained 
using BFS and HK76 algorithms are identical within er- 
ror bar. 

The mean degree of the system anisotropy has been 
calculated as 



2=1 



(9) 



where Si = (R% — Rf)/Ri. Here, R% and Rf are radii of 
gyration of cluster i in y and x directions, respectively, 
Ri is its mean radius of gyration, N c is a total number 
of clusters, Ni is a number of filled sites in the cluster z, 
and N t is a total number of the filled sites. 



III. RESULTS AND DISCUSSION 
A. Non-universality of intersection points R* 

The value of the percolation probability or percolation 
cumulant at the intersection point R* may be impor- 
tant characteristic representing the universality class (l3l |. 
Figure 0] presents examples percolation probability Rl 
versus /c-mers concentration p for isotropic systems, s = 
0, and different values of k and L. 

For isotropic problem the position of intersection point 
remained unchanged within precision of estimation, be- 
ing R* ~ 0.90 for all k within the interval between 1 
and 512. This behavior is rather similar to that observed 
for percolation problem of /c-mers with open boundary 
condition [8j. For the criterium and, the same values 
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Figure 4. (Color online) Probability curves for isotropic sys- 
tems, s — 0, and different values of k and L. Arrows indicate 
the intersection points. 



of R* ~ 0.3 have been observed for the different length 
of k-mers ranging between k = 1 and k = 25. Thus, 
universality of intersection points R* has been observed 
for the systems with different boundary conditions (pe- 
riodical and open) and it may indicate the conserving of 
universality class irrespective of the size of /c-mers. 

However, such universality of intersection points R* 
has been not observed for anisotropic systems. Figure [5] 
presents examples of percolation probability R versus k- 
mer concentration p for k = 32 and different values of s 
and L. At fixed value of k the position of intersection 
point i?* continuously decreased with increasing of s. 

The more detailed studies have shown that for 
anisotropic systems the position of intersection points 
R* also dependents on the value of k (Figure [6]). For 
the completely ordered systems, s = 1, the value of R* 
decreased monotonically and became close to for larger 
sizes of /c-mers. This observation may reflect the contin- 
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Figure 5. (Color online) Probability curves for k = 32, s 
0.0,0.7, 1.0. Arrows indicate the intersection points. 



uous change of universality class and correspondences to 
previously reported data for the completely ordered sys- 
tems with open boundary conditions [13|- For partially 
ordered systems, the similar effect of /c-mers length on 
the value of R* has been observed (Figure [6]). 
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Figure 6. (Color online) Intersection point of percolation 
probability R* versus k at s = 0.0, 0.7, 1.0. 

Thus, orientation of k-mers affected the universality 
class of this percolation problem and it has been con- 
served only for the isotropic systems (s = 0), where uni- 
versality is the same for the different length of /c-mers. 
It can be speculated that this violation of universality 
can reflect the effect of the system anisotropy. This 
anisotropy has been maximally denominated for the com- 
pletely ordered systems (5 = 1) where the effect of the k- 
mer length on the value of R* is maximal (Figure [6]). The 
more detailed analysis shown that the structure of perco- 
lation clusters is strongly depends upon fc, they have been 



elongated along vertical direction and degree of elonga- 
tion increased as length of /c-mers increased (Figure [TJ. 
Moreover, the mean degree of system anisotropy 5 calcu- 
lated using Eq. [9] is dependent on k-mei concentration p, 
length k and order parameter s. 

Figure [Hfci presents examples of S versus order parame- 
ter s at different fixed concentrations p and fixed length 
of /c-mer, k = 32. The size of lattice is relatively large, 
L = 4096, so, the finite size effects are rather small. For 
isotropic systems (at s = 0), the value of 5 is always 
zero and it is maximal for completely ordered systems 
(at s = 1). At small values of p the relation between S 
and s is nearly linear. With increasing of p and fixed s 
the value of S decreased, however, it noticeably dropped 
above percolation threshold and became practically zero 
in the vicinity of jamming concentration. E.g., the con- 
centration of p = 0.50 is above the percolation threshold 
for the systems with order parameter s below ~ 0.8 and 
here, the S(s) dependence noticeably deviates from near 
linear (Figure [Ht)- Figure [SJd presents examples of S ver- 
sus order parameter s at different fixed length of /c-mer, 
k = 32 and the concentrations that has corresponded 
to the percolation transitions for the given systems. The 
value of fc-mer length k strongly affected the mean degree 
of the system anisotropy S at the percolation transition. 
E.g., for dimers, k — 2, the value of 5 is rather small 
in the whole range of s between and 1, however, with 
increasing of A:, the 5(s) became more noticeable and we 
believe that they can transfer into the near-linear of type 
S ~ s in the limit of large k-mer length, k — » oo. 



B. Dependence of percolation threshold p c versus 
order parameter s 

The p c (s) dependencies for /c-mers of different length 
(fc = 2 . . . 128 are presented in Figure [9] For complete- 
ness, the precise numerical information is also collected 
in Table |U In addition the Table [TT] presents rougher es- 
timations for k = 256 and k = 512 for isotropic (s = 0) 
and completely ordered (s = 1) systems. 

The obtained data evidence that the increase of sys- 
tem ordering always results in increase of p c value. Such 
behavior correlates with theoretical data obtained in the 
systems of partially oriented penetrable rods (37l-[40j and 
experimentally studied effect of carbon nanotube align- 
ment on percolation in polymer composites [41]. 

Figure [10] presents examples of p c versus k dependen- 
cies obtained for different values of s in this work (1), as 
well as data presented earlier for the isotropic (s = 0) 
and completely ordered (s = 1): (2) @, (3) 0, (4) (HJ, 
(5) H and (6)13. 

For completely ordered systems, i.e., at s = 1, the per- 
colation threshold p c monotonically decreased as value of 
k increased. Recently, the similar behavior for k from 1 to 
12 with the asymptotic limit of = p c (k — > oo) ~ 0.54 
has been reported [13] . 

The analysis have shown that the data obtained in our 




Figure 7. Examples of wrapping clusters incipient in vertical direction for completely ordered system (s = 1.0) for different 
length of /c-mers. The size of a square lattice is L = 128k. Periodical boundary conditions. 
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Figure 8. Mean degree of the system anisotropy S versus order parameter s: a) at different concentrations p and the fixed 
length of a /c-mer, k = 32; b) at different k-mei length k and concentrations that have been corresponded to the percolation 
transition for given systems. The size of lattice is L = 4096 and the data have been averaged on 100 independent runs. 



Table I. Percolation threshold p c versus order parameter s for /c-mers of different length k. 



s 


k = 2 


k = 4 


k = 8 


k = 16 


k = 32 


k = 64 


k = 128 


0.0 


0.5619 


0.5050 


0.4697 


0.4638 


0.4748 


0.4928 


0.5115 


0.1 


0.5621 


0.5056 


0.4702 


0.4644 


0.4751 


0.4930 




0.2 


0.5627 


0.5067 


0.4717 


0.4656 


0.4763 


0.4936 




0.3 


0.5638 


0.5090 


0.4742 


0.4677 


0.4777 


0.4948 




0.4 


0.5653 


0.5124 


0.4777 


0.4708 


0.4802 


0.4964 




0.5 


0.5672 


0.5167 


0.4825 


0.4751 


0.4834 


0.4993 




0.6 


0.5698 


0.5224 


0.4890 


0.4807 


0.4879 


0.5025 




0.7 


0.5728 


0.5296 


0.4977 


0.4883 


0.4939 


0.5074 




0.8 


0.5765 


0.5389 


0.5092 


0.4987 


0.5021 


0.5132 




0.9 


0.5809 


0.5510 


0.5251 


0.5140 


0.5142 


0.5210 




1.0 


0.5862 


0.5672 


0.5526 


0.5442 


0.5397 


0.5376 


0.5366 
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Figure 9. (Color online) Percolation threshold p c versus order 
parameter s for /c-mers of different length. 



Table II. Estimations of percolation threshold p c for /c-mers 
of ladge length k. 



s 


k 256 


k 512 


0.0 
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0.5485 


1.0 


0.535 
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Figure 10. (Color online) Percolation threshold p c versus k- 
mer length k at different values of order parameter s. Here, 
the different data are presented that have been obtained in: 
(1) this work, (2) H, (3) 0, (4) E2, (5) @ and (6) [H]. The 
dashed lines have been obtained by least square fitting of the 
data points using the Eqs. 1101111 



work may be rather well fitted by the power function 

p c = a 1 /k^+p™, (10) 

where p* = 0.533±0.001, a x = 0.088±0.003, a ± = 0.72± 
0.04 and r 2 = 0.998 for the coefficient of determination. 
Note, that for completely ordered penetrating 



anisotropic objects and continuous problem the excluded 
volume theory predicts the absence of noticeable de- 
pendence of the percolation threshold on aspect ratio 
k |37|, l38| . In our lattice problem, observed effect of 
p c (k) dependence may reflect the influence of the lattice 
discreteness on the percolation threshold. 

In our problem at s = 1, the formation of percolation 
cluster reflects the mode of connectivity between verti- 
cally oriented one dimensional chains of k-meis. It may 
be assumed that in the limit of k — >• oo the connectivity 
of two /c-mers in the neighbor vertical lines at their end 
sites is sufficient for a formation of percolation cluster 
with minimal concentration of p = 0.5, that is close to 
the numerically estimated value of = 0.533 ±0.001. 

In contrast, for partially ordered systems, i.e., at s < 1, 
the percolation threshold p c is a nonmonotonic function 
of k and for a certain length of k-meis k = k m a minimum 
of p c has been observed (Figure [TO]) . In total, the data 
obtained in this work for isotropic systems (i.e., at s = 0) 
have been in good correspondence with previously pub- 
lished data 0, [1, EH, [Hi with only exception to those 
obtained for the very long /c-mers (k > 64) in |42|. This 
inequality may reflect the relatively moderate size of lat- 
tices that has been used in [42] (L < 2500) whereas in our 
simulations the maximum size of the lattice L ~ 100 A:, as 
a rule. In any case, our data confirmed the conclusion by 
Leroyer and Pommiers [5] and Kondrat and P§kalski [6] 
about the presence of minimum at the p c versus k de- 
pendence. For the disordered systems the position of the 
minimum, fc m , is dependent on the value of s, e.g., it 
is k rn ~ 13 at s = 0, k rn ~ 16 at s = 0.7, k rn ~ 22 
at s = 0.9, and it seems that k m —> oo in the limit of 
s —> 1 (Figure [10]). Note, that the asymptotic limit of 
Pc = Pc(k — » oo) ~ 0.461 derived in [8] for s = in 
fact is very close to the value of p c at point of minimum, 
k rn — 13. 

It is attractive to speculate that extremal p c versus be- 
havior for partially ordered systems may reflect the com- 
petition of the two different effects influencing the value 
of percolation threshold. We tried to fit the obtained 
data for isotropic system (s = 0) using the function 

p c = a /k ao + Mog 10 k + c (11) 

and had obtained the following numerical estimations for 
the parameters ao = 0.36 ± 0.02, ao = 0.81 ± 0.12, b — 
0.08 ± 0.01, c = 0.33 ± 0.02, and r 2 = 0.991 for the 
coefficient of determination. 

It is remarkable, that exponents ao = 0.81 ±0.12 and 
ql\ = 0.72 ± 0.04 are practically the same for s = and 
5 = 1, respectively, and it may reflect the same effect 
of the discreteness on the percolation at the relatively 
small k (< 10). For disordered systems, the logarithmic 
increase of p c at large values of k (Eq fTT]) may reflect the 
tendency of k-mers for stacking, or formation of square- 
like blocks, especially at large values of k. Such blocks 
of vertically and horizontally oriented /c-mers are typi- 
cal forpartially ordered systems in jamming configura- 
tions [12j], however, they are also important at the perco- 



(a) (b) (c) 

Figure 11. (Color online) Examples of percolation configurations of /c-mers (k — 128) on a square lattice of size L — 4096 
at s = 0. Vertical and horizontal orientations are represented by different gray levels in printed version and in red and blue 
in online version; empty sites are labeled black and sites of percolation cluster are labeled white. Here, (a) shows the whole 
lattice; (b) shows the magnification of the pattern (a) in the central square; (c) shows the magnification of the pattern (b) in 
the central square. 



lation threshold. Examples of k-mer patterns (k = 128) 
in the percolation point are presented in Figure [11] for 
isotropic system (s = 0). The sequential magnification 
of the system has shown the presence of rather com- 
pact blocks of vertically and horizontally oriented k-mers 
that have been connected into the percolating structure 
by overhanging of /c-mers. The numerical studies have 
shown that for the ideal blocks, i.e. k x k squares the per- 
colation concentration increased and jamming concentra- 
tion decreased as k value increased [43] and above certain 
critical value of k no percolation has been observed. In 
this situation even at the saturation coverage (jamming) 
where no more object can be placed without any overlap 
there exist only finite clusters of k x k squares. We can 
assume the similar mechanism that governs the observed 
p c oc k increasing of percolation threshold. 

We checked the validity of conjecture of Vandewalle et 
al. [3] about the constancy of the ratio of percolation and 
jamming concentration p c /Pj for disordered (s = 0) and 
completely ordered 5 = 1 systems (Figure H2|) . The val- 
ues of jamming concentration pj have been taken from 
our previously published work [12] . For completely or- 
dered systems (s = 1) this ratio initially increased and 
became practically constant, p c /Pj — 0.715, at relatively 
large length of fc-mers, k > 8. 

For isotropic systems (s = 0) this ratio is approxi- 
mately constant only small values of k (k = 2 ... 8) and 
for larger /c, k = 16 . . . 256 the p c /Pj increased propor- 
tionally to log 10 k: 

Pc/Pj = frlogio k + C i ( 12 ) 

where b = 0.119 ± 0.003 and c = 0.513 ± 0.006 are the 
constants. 

Thus, the constancy of the ratio p c /Pj fulfills for the 
completely ordered systems (s = 1). It may reflect the 
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Figure 12. (Color online) Ratio of percolation and jamming 
concentration p c /pj versus k-mei length k for disordered (s = 
0) and completely ordered s = 1 systems. The dashed lines 
for s = has been obtained by least square fitting of the data 
points using the Eqs. [12] 



similar influence of the discreteness of the lattice on the 
both jamming and percolation. On the other hand, the 
non-constancy of the ratio p c /Pj for isotropic systems 
(s = 0) may reflect the different influence of the stack- 
ing on the jamming and percolation. For this case the 
approximation of the data presented in Figure [12] gives 
p c /pj ~ 1 at fc ~ 1.2 x 10 4 . So, we can suppose that 
for a very long /c-mers the percolation may be lost in 
close analogy with similar behavior observed for k x k 
squares |43|. 
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IV. CONCLUSION 

In this paper, the percolation behavior of partially or- 
dered linear fc-mers on torus (square lattice with periodic 
boundary conditions) has been investigated by computer 
simulations. The length of a k-mer varies from 1 to 512 
and different lattice sizes up to L = 1024/c are used. The 
relaxation random sequential adsorption model [12] has 
been used to place the k-mers on a lattice. The alignment 
degree is characterized by order parameter s = ... 1 : 
s = for isotropic system and s = 1 for perfectly aligned 
system. The behavior of percolation cumulant at the in- 
tersection point R* has been studied in details in depen- 
dence on fc, s and L. For isotropic problem the value of 
position of intersection point remained unchanged within 
precision of estimation, being R* ~ 0.90 for all studied 
length of k-mer s. The universality of intersection points 
R* (i.e., independence of R* on k) has been observed 
only for isotropic systems, s = 0. This universality sug- 
gests that R* can be derived from the work (3]J not only 
for topological percolation but also for physical one. For 
anisotropic systems this universality is violated and the 
value of R* is dependent upon k and s. One can suppose 
that this violation can reflect the effect of the system 
anisotropy. 

The increase of system ordering always results in in- 
crease of percolation threshold p c . The dependencies 
of p c (k) for completely ordered (5 = 1) and partially 
ordered (s < 1) systems are obviously different. For 



completely ordered systems the percolation threshold p c 
monotonically decreased as k increased. The power law 
relation p c oc l/k ai (a x = 0.72 ± 0.04) probably reflects 
effects of the lattice discreteness. For partially ordered 
systems the percolation threshold p c is always a non- 
monotonic function of k and for a certain length of k- 
mers k = fc m a minimum of p c has been observed. It has 
been assumed that this behavior may reflect the com- 
petition of the lattice discreteness (that is dominant at 
small values of k) and the tendency of /c-mers for stack- 
ing, or formation of squarelike blocks (that is dominant 
at large values of k). For completely ordered systems 
(s = 1) the ratio of percolation and jamming concen- 
tration Pc/Pj is practically constant (p c /Pj — 0.715, at 
k > 8). This behavior evidently reflects the presence of 
some universal connection in the geometry of percolation 
and jamming Q. For isotropic systems (s = 0) this ratio 
is not constant and increased proportionally to log 10 k. 

Our simulations suggest that for s = the percolation 
may be lost at k }t 1.2 x 10 4 . Additional investigation 
of percolation with extreme long objects should be per- 
formed in future to confirm or reject this prediction. 
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